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ABSTRACT 


The  work  done  for  the  project  is  reported  here  in  three  parts.  The 
first  is  the  analysis  for  minimization  of  the  temperature  distortion  due  to 
the  thermocouple  cavity.  The  error  is  minimized  or  reduced  to  zero  by  opti- 
mizing a combination  ot  cavity  diameter  and  depth  and  the  thermocouple 
material  and  size.  The  second  is  the  refinement  of  the  presently  available 
computer  programs  for  prediction  of  the  surface  temperature  and  heat  flux 
at  the  inner  surface  of  the  pipe  by  inverting  the  temperature  measured  by 
an  interior  probe  close  to  the  heated  surface.  The  refinement  is  achieved 
by  using  the  double  precision  format  in  the  program  and  adapting  the  dimen- 
sionless formulation.  The  third  is  to  study  the  inversion  solution  for  a 
large  time  duration  of  a time  dependent  surface  heat  flux.  The  solution  is 
obtained  by  the  method  of  Laplace  transform  and  the  convolution  integral. 

Each  of  the  above  three  subjects  is  reported  as  a part  of  the  present  report. 
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FAFT  r ON  MINIMIZATION  OF  TEMPERATURE  DISTORTION 
IN  THE  THERMOCOUPLE  CAVTTY 

I.  INTRODUCTION 

A direct  measurement  of  transient  surface  temperature  and  heat  flux 
is  often  difficult.  For  example,  a surface  involves  two  modes  of  heat 
transfer,  say,  radiative  and  convective  heat  transfer.  In  this  case  If 
the  measuring  probe  has  a different  radiative  property  than  that  of  the 
surface,  erroneous  measurements  will  result.  A piston  or  projectile 
sliding  over  the  cylindrical  surface  is  another  case  where  the  direct 
measurement  at  the  surface  is  difficult.  A surface  involving  melting  or 
ablation  is  also  difficult  to  make  direct  measurement.  Therefore,  in- 
direct estimation  by  inverting  the  temperature  history  inside  the  heat 
conducting  solid  as  measured  by  a thermocouple  is  often  used  for  pre- 
diction of  the  surface  temperature  and  heat  flux.  Beck  [1],  Hernning 
and  Parker  [2],  Frank  [3],  Imber  and  Khan  [4],  Stolz  [5],  Chen  and 
Thomsen  [6]  have  developed  different  inversion  solutions  for  this  pur- 
pose. All  of  these  solutions  assumed  that  the  cavity  drilled  into  the 
solid  does  not  distort  the  true  temperature  distribution.  Thus,  it  is 
important  that  the  temperature  measurement  by  an  interior  probe  is 
accurate  and  involves  least  dirtortion  or  error.  Theoretically  Beck  [7], 
Masters  and  Stein  [8],  Burnett  [9],  and  Chen  and  Li  [10]  studied  the  distor- 
tion of  the.  temperature  field  in  the  present  of  thermocouple  and  its 
cavity.  Experimentally  Chen  and  Danh  [11]  showed  that  appreciable  dis- 
tortion, say  10%,  of  temperature  field  may  exist  for  a normal  implant 
of  the  thermocouple  into  s solid  body.  From  studies  of  Chen  and  Li  [10] 
and  Beck  [7]  they  found  that  with  a proper  combination  of  the  thermocouple 
cavity  diameter,  its  depth,  and  the  thermocouple  material  and  its  diameter 


the  distortion  of  temperature  field  with  respect  to  space  or  time  can 
be  minimized  if  not  eliminated.  In  this  report  we  study  the  optimum 
combination  of  these  parameters  such  that  at  a given  situation  one 
knows  what  is  the  best  combination  to  use  and  what  is  the  magnitude  of 
the  temperature  distortion. 


II.  FORMULATION  OF  PROBLEM 


In  the  present  study  we  consider  a disk  depicted  in  Figure  1 which 
has  a thickness  D and  is  drilled  a cavity  of  a diameter  to  a depth 
of  e distance  from  the  heating  surface.  The  heat  flux  Q is  assumed  to 
be  constant  and  the  upper  surface  of  the  disk  is  assumed  to  be  insulated. 

A thermocouple  of  a diameter  d^  is  then  welded  on  the  cavity  base. 

Furthermore,  the  disk  may  be  thought  to  approximate  a cut  out  from  a 
hollow  cylinder  if  the  radius  of  the  disk  is  small  compared  with  the  radius 
of  the  cylinder.  The  diameter  of  the  disk  is  choosen  to  be  2D  outside  which  the 
temperature  distortion  due  to  the  thermocouple  cavity  becomes  negligiable. 

For  this  to  be  true  one  needs  to  restrict  the  ratio  of  cavity  diameter 
to  the  disk  diameter  d/2D  be  small.  The  unfilled  cavity  can  be  air  or 
insulating  material. 

The  basic  idea  to  minimize  or  to  eliminate  the  temperature  distortion 
is  based  on  a proper  choice  of  the  thermocouple  size  and  the  material 
which  has  a higher  thermal  conductivity  than  that  of  the  disk  so  as  to 
conduct  more  heat  away  at  the  cavity  base  balancing  the  insulation  effect 
of  the  insulator  in  the  cavity. 

Let  X and  Y be  respectively  the  coordinate  along  and  normal  to  the 

i 

heated  disk  surface  and  the  Y axis  coincide  vrith  the  axis  of  the  cavity. 
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The  thermal  conductivity  is  assumed  to  be  constant  and  the  temperature 
distribution  is  axisymmetrical . The  governing  equations  for  the  tran- 
sient heat  conduction  in  dimensionless  form  are: 
for  the  disk  (subscript  "1",  see  Figure  2) 


30  3_0  1 30  3 0 

9t  a 2 + x 3x  + . 2 

3x  3y 


for  the  insulating  material  in  the  cavity  (subscript  "2") 


3M9  _ _2 
3x  a. 
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for  the  thermocouple  (subscript  "3") 
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where  X = a^t/D  is  the  dimensionless  time,  x = X/D  the  dimensionless 
radial  coordinate,  y = Y/D,  the  dimensionless  distance  normal  to  the 
heated  surface,  a^,  a^,  and  are  respectively  the  thermal  diffusivity 
oi  i lies  disk,  insulating  material,  and  the  thermocouple.  The  dimension- 
less temperature  0 is  defined  as  Tk^/QD  where  T is  the  temperature  above 
the  initial,  uniform  temperature. 

The  intial  temperature  of  the  disk  is  then 


0 (x,  y,  0)  = 0 


The  boundary  conditions,  see  Figure  2,  are: 


a constant  heat  flux  at  the  lower  surface. 


y = 0 


1 


the  Insulation  at  the  upper  surface. 


y = 0 


0 


). 


(5) 
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’-he  zero  temperature  distortion  at  the  edge  of  the  disk, 


x = 1 
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and  the  axisymmetric  condition  at  the  cavity  axis 
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= 0 


39 

3x 


= 0 


x=0 


(8) 


The  condition  (7)  of  the  zero  temperature  distortion  was  verified  by 
Chen  and  Li  [10]  in  their  early  calculation  when  the  cavity  diameter  is 
one  tenth  of  the  disk  diameter.  In  addition  to  the  above  boundary 
condition  the  temperature  and  heat  flux  at  the  interface  of  the  disk, 
thermocouple,  and  insulating  material  are  taken  to  be  continuous. 

Ill  ANALYSIS 

There  are  five  parameters  that  can  be  varied  for  the  present  anal- 


ysis. They  are  (a)  the  dimensionless  distance  from  the  base  of  the 
cavity  to  the  heated  surface  e/D,  (b)  the  size  of  the  cavity  d/D, 

(c)  the  ratio  of  the  thermocouple  diameter  to  that  of  the  caviLy,  dt/d. 


(d)  the  thermal  conductivity  ratio  and  which  comes  from  the 

continuity  of  heat  flux  at  interfaces.  (e)  the  ratio  of  the  product  of 
density  and  specific  heat  Pj^/p]^  (or  equivalently  to  the  ratio  of 
of  thermal  diffisivity  or  o^/c^).  The  subscripts  1,  2 and  3 denote 

the  disk,  insulation  and  thermocouple  materials. 

Because  of  the  complexity  of  the  geometry  and  the  multicity  of 
material  the  method  of  finite  element  technique  as  discus; ed  by  Wilson 
and  Nikel  [12]  is  adapted  with  the  aid  of  a computer  program  developed 
by  Wilson  [13].  The  present  problem  is  subdivided  into  finite  element 
as  required  by  the  method,  see  Figure  2.  The  dimensionless  pie  section 
is  subdivided  into  121  finite  elements  with  12  dividing  lines  on  both 
coordinates.  Each  element  is  defined  by  four  nodal  points  where  nodal 
points  are  denoted  by  intersections  of  the  dividing  lines  and  numbered 
as  shown  in  Figure  2.  The  material  property  corresponding  to  each 
element  is  then  assigned  to  the  program  developed  by  Wilson  [13] 

The  solution  at  each  nodal  with  respect  to  time  is  then  obtained. 

For  numerical  calculation  three  typical  values  of  the  distance 
from  heated  surface  to  the  base  of  the  cavity  e/D  are  chosen  to  be  0.04, 
0.1  and  0.2.  The  cavity  diameter  is  fixed  at  one  tenth  of  the  disk 
diameter.  The  thermocouple  to  cavity  diameter  ratfo  d^/d  is  made  to 
vary  0,  0.2,  0.4,  0.6,  0.8  and  1.0.  Regarding  the  range  of  the  ratio 
of  the  thermal  conductivity  and  the  ratio  of  the  product  of  density  and 
specific  heat  we  surveyed  these  ratios  for  the  commonly  used  thermocouples 
in  Table  1 [14]  and  plotted  in  Figure  3.  Figure  3b  shows  the  ratio  of 


thermocouple  conductivity  to  that  of  the  disk  versus  the  ratio  of 

density  - specific  heat  product  Pj^/P-j^j  where  the  value  of  the  conduc- 
tivity is  taken  to  be  the  average  value  between  200  and  800°K.  One  secs 
that  for  most  practical  situations  the  ratio  of  density  - specific  heat 
product  P3C3^P^C^  *s  approximately  constant  at  0.7  except  when  the  con- 
ductivity ratio  k ^ /k ^ is  small.  Thus  the  value  of  p3C3/P-LC3  and 
for  calculation  are  chosen,  as  shorn  in  triangular  symbols  of  Figure  3b, 
to  cover  the  practical  range.  The  corresponding  value  for  ^2C2^1C"'  anc^ 
K2^K1  ^°r  *nsu-*-ati°n  material  are  chosen  to  be  fixed  at  0.5  and  0.005 
which  is  a typical  value  for  Teflon  insulating  material  and  is  also  ap- 
proximately the  order  of  magnitude  for  air. 

IV.  RESULTS  AND  DISCUSSIONS 

Numerical  results  of  the  calculations  are  presented  in  Tables  2 to 
4 and  Figures  4 to  9.  The  percentage  error  of  temperature  is  defined 
as  the  distorted  temperature  divided  by  a reference  temperature  defined 
by  QD/k^.  Tables  2 to  4 give  the  percentage  error  of  temperature  dis- 
tortion) as  a function  of  time  for  different  values  of  the  parameters 
djd  (0  to  1.0),  (0.5  to  10),  P3C3^P^C3  (0.5  to  1.8)  and  e/D 

(0.02  to  0.1). 

Figures  4,  5 and  6 show  the  three  typical  temperature  distributions 
in  the  steel  disk  near  the  thermocouple  junction  for  the  case  e/D  = 0.06 
d/2D  = 0.1  at  the  time  r = 0.08.  Figure  4 is  the  temperature  distribu- 
tion when  the  cavity  is  filled  entirely  with  the  insulation  material 
(Teflon  "3"  k.,/^  - 0.005,  ~ 0*5).  The  dimensionless  isotherms 
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Tic  /QD  shows  the  distortion  of  temperature  distribution.  One  sees  that 
1 

the  insulation  effect  on  the  heat  transfer  near  the  cavity  base  not  only 
creates  a much  higher  junction  temperature  of  Tk^/QD  = 0.342  than  the 
undistorted  one  of  0.255  at  the  edge  of  the  disk  giving  an  error  of  8.7% 
but  also  causes  a hot  spot  at  the  heating  surface  with  a higher  tempera- 
ture  of  Tk^/QD  = 0.374  over  the  undistorted  one  of  0.311.  On  the  other 
hand  the  temperature  distribution  in  the  insulation  material  is  much 
lower  than  the  true  temperature  creating  a large  temperature  gradient 
at  the  base  of  the  cavity.  Figure  5,  contrary  to  Figure  4,  is  the  tem- 
perature distribution  when  the  cavity  is  completely  filled  with  thermo- 
couple material  whose  thermal  conductivity  is  ten  times  larger  than  the 
disk  material  e.g. , copper  versus  steel.  Now  over-conduction  of  heat 
by  the  thermocouple  has  created  a cold  spot  at  the  base  of  the  cavity 
giving,  Tk^/QD  of  0.17  versus  the  undistorted  one  of  0.255  with  an 
error  of  8.5%  as  well  as  at  the  heating  surface  with  Tk^/QD  of  0.236 
versus  0.311.  The  temperature  distribution  in  the  thermocouple  now 
becomes  higher  than  the  undistorted  one  in  the  disk.  By  properly 
choosing  the  ratio  of  thermocouple  diameter  to  that  of  the  cavity  one 
may  minimize  these  distortions  of  the  temperature  response  at  the  base 
of  the  cavity.  This  is  shown  in  Figure  6 where  the  thermocouple 
diameter  d^  is  choosen  to  be  0.4  of  the  cavity  diameter  d with  = 10* 

e.g.  copper-steel  combination.  Figure  6 shows  that  distortions  at  the 
base  of  the  cavity  and  at  the  heating  surface  are  almost  eliminated  giving 
the  temperature  Tk^/QD  of  0.245  and  0.320  with  respectively  the  error  of 
1%  and  0.9%. 
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In  order  to  examine  the  details  of  the  distorted  temperature  response 
at  the  base  of  the  cavity  we  tabulated  the  results  in  Tables  2,  3 and  4 
and  plotted  the  error  percentage  as  function  of  the  ratio  of  the  thermo- 
couple diameter  to  that  of  the  cavity  for  different  cavity  depth,  time, 
and  thermal  conductivity  in  Figures  7,  8 and  9.  In  these  Figures  the  pc 
ratio  ranges  from  0.7  to  1.3  which  covers  most  of  the  practical  applica- 
tion. 

In  the  case  when  the  pc  ratio  is  equal  to  or  less  than  one  the  error 
in  temperature  response  at  the  base  of  the  cavity  in  these  figures  are 
all  positive  or  overheat  for  k^/k^  < 1 • This  is  because  the  thermal 
conductivity  of  the  thermocouple  is  less  than  that  of  the  disk  and  the 
heat  capacity  pc  of  the  thermocouple  is  also  small.  Therefore,  no  extra 
conduction  of  heat  can  be  achieved  by  the  thermocouple  to  compensate  for 
the  blocking  of  the  heat  transfer  by  the  insulation  material  in  the 
cavity.  On  the  other  hand,  if  > 1 the  error  of  temperature  varies 

from  positive  value  for  d^/d  = 0 to  some  negative  value  as  d^/d  approaches 
1.  Thus  for  k^/k^  > 1 a properly  chosen  combination  of  thermocouple 
and  insulation  material  can  minimize  the  error.  For  example,  in  Figure  7 
a combination  of  k^/k^  = 10,  P3C3^pici  = 0-75,  e/D  = 0.1  and  d^/d  = 0.5 
produces  almost  negligible  error.  This  combination  shown  to  be  optimum 
at  T = 1.0  in  Figure  7 is  also  optimum  for  other  time  periods  (see  Table 
2).  Therefore  once  an  optimum  combination  of  parameters  is  chosen  it  is 
valid  throughout  the  entire  transient  period  of  an  experiment.  From 
Figures  7,  8 and  9 one  can  also  see  that  the  optimum  ratio  of  djd  which 
gives  zero  temperature  error  decreases  as  the  ratio  increases. 


i 


i 


This  implies  that  for  the  thermocouple  with  a larger  thermal  conductivity 
a smaller  diameter  is  sufficient  to  eliminate  the  temperature  distortion. 
The  result  shown  in  Figures  7,  8 and  9 can  in  general  be  adapted  for  use 
in  practical  application  to  choose  the  size  of  thermocouple  and  cavity, 
the  thermocouple  material  and  the  depth  of  the  cavity  to  be  drilled. 

As  mentioned  earlier  that  when  < 1 and  p3c3^pici  - * the  error 

of  the  temperature  response  at  the  base  of  the  cavity  are  all  positive. 
However  we  found  (see  Tables  2.2,  3.2  and  4.2)  that  if  Pj^/Pj^j  ratio 
is  made  large  enough  during  the  transient  period  the  error  of  temperature 
response  at  the  base  of  the  cavity  may  indeed  become  negative  even  when 
k^/k^  < 1.  Physically  although  the  thermal  conductivity  of  the  thermo- 
couple is  smaller  than  that  of  the  disk  material,  but  with  a larger 
heat  capicitance  *>2Z3^1C1  > * the  thermocouple  is  still  capable  of 
absorbing  extra  heat  flux  and  hence  eliminates  the  temperature  distor- 
tion at  the  cavity  base  during  the  transient  period.  To  illustrate  this 
fact  we  examine  Figure  7 (or  see  Table  4.2)  for  the  data  of  K^Ki  = 1 
and  p^c^/p^c^  = 1*3.  One  sees  that  when  d^/d  = 1 the  temperature  dis- 
tortion can  indeed  be  negative.  Therefore,  if  d^/d  are  choosen  between 
0.8  and  1 the  error  can  be  minimized.  However  one  must  keep  in  mind 
that  the  elimination  of  error  by  heat  capacitance  can  work  during  the 
transient  period  only  for  once  a steady  state  conduction  is  established 
the  heat  capacity  pc  will  no  longer  have  any  effect  and  over  heating 
at  the  cavity  base  eventually  will  develop.  This  can  be  seen  best  from 
the  governing  equation  (1)  that  at  steady  state  the  unsteady  term  which 
contains  pc  product  is  zero  and  is  not  a parameter  affecting  the  dis- 


tortion. 


Another  important  fact  that  should  be  mentioned  is  that  in  general 
the  optimum  choice  of  df/d  ratio  for  given  k^/k^  and  Pc  ratio  does  not 

vary  very  much  with  the  variation  of  e/D  ratio.  The  insensitivity  of 

( 

the  optimum  d£/d  ratio  to  the  c/D  ratio  ranging  from  0.02  to  0.1  means 
that  the  distoi cion  of  the  temperature  is  insensitive  to  the  cavity  depth 
or  the  thickness  of  the  disk.  This  fact  was  already  pointed  out  by  Chen 
and  Danh  [11]  in  their  experiment  that  the  temperature  distortion  at  the 
base  of  the  cavity  is  more  sensitive  to  the  variation  of  the  cavity 
diameter  than  the  depth  of  the  cavity  drilled. 

As  an  example  of  a practical  applicacion,  let  us  consider  a measure- 
ment of  the  transient  temperature  response  of  an  engine  block  made  of 
aluminum.  From  Figure  3b  we  know  that  aluminum  has  high  thermal  conduc- 

I 

tivity.  Therefore  copper-constantan  thermalcouple  which  has  a higher 
thermal  conductivity  than  aluminum  should  be  chosen.  For  this  material 
combination  we  have  kt^/k^  = 1*69  p^c^/p^c^  = 1. 3.  Now  if  the  thermocouple 

cavity  is  drilled  such  that  e/D  =0.1  then  from  Figure  7 interpolating 
between  cc^/ic^  = 2 and  1 for  P3c3/p^c^  = 1-3  we  find  that  the  optimum 
df/d  for  ~ 1-69  is  approximately  0.7 

One  disadvantage  of  invoking  finite  element  analysis  is  that  the 
result  does  not  give  a clear  functional  relation  among  the  parameters 
involved.  In  an  attempt  to  obtain  a simple  and  useful  relation  to 
relate  the  various  parameters  we  note  the  following  fact  and  result: 

(a)  the  optimum  d /d  ratio  for  zero  temperature  distortion  is  a strong 
function  of  and  pc  ratio  but  is  relatively  insensitive  to  the  i/D 

ratio,  (b)  from  the  theoretical  reasoning  the  d^/d  ratio  is  independent 
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of  pc  ratio  if  the  problem  is  steady  state.  A simple  steady  one  dimen- 
sional analysis  in  which  the  thermocouple  and  thi  : nsulati  ri  material  in 

the  cavity  is  made  to  conduct  the  same  amount  of  heat  that  would  be 
transferred  without  the  cavity  gives  the  relation 

d^/d  = ✓(*  - ic2^^k3  - 

Using  the  above  equation  as  a base  we  find  that  for  the  transient  heat 
conduction  as  calculated  by  the  finite  element  method  the  following 
equation  (10)  correlates  very  well  with  the  optimum  dt/d  ratio. 

d (_ / d = (pjC.j/p^c^)0'  3 / (k^  - K.p/ - <2)  (10) 

Equation  (10)  gives  an  error  or  distortion  of  no  more  than  two  per- 
centage r ->ints.  In  practice  equation  (10)  may  be  used  as  a rule  of 

thumb. 

V.  CONCLUSION 

An  analysis  of  the  temperature  distortion  caused  by  the  cavity 
drilled  into  a disk  to  accommodate  the  thermocouple  has  been  studied. 

The  calculation  is  carried  out  for  the  case  of  constant  heat  flux.  It 

is  shown  that  the  temperature  at  the  base  of  the  cavity  distorted  from 
that  without  a cavity  can  be  eliminated  by  a properly  chosen  combina- 
tion of  the  ratio  of  the  thermocouple  diameter  to  the  cavity  diameter, 
d^/d  and  the  thermocouple  material  k^/k^.  The  optimum  ratio  of  d^/d 

can  be  found  from  Figures  7,  8 and  9 or  Tables  2,  3 and  4,  or  approxi- 

mately from  equation  (10).  As  a rule  the  thermocouple  must  be  chosen 


to  have  a higher  thermal  conductivity  than  that  of  the  heat  conducting 
solid.  The  cavity  diameter  should  be  as  small  as  practically  possible. 
For  the  case  of  time  dependent  surface  heat  flux  the  present  result  may 


be  also  used  as  a general  guide. 
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THERMAL  CONDUCTIVITY 


INSULATING  MATERIAL  ( TEFLON 


rigure  4 Temperature  Distribution  with  Insulation  Material  Filled 


ALL  THERMOCOUPLE  MATERIAL 


Figure  5 Temperature  Distribution  with  Thermocouple  Material  Filled 


THERMOCOUPLE  ( 3 ) 

/—  INSULATION  MATERIAL 


Temperature  Distribu^-on  with  Thermocouple  Material  Partially 


ERROR 


Ratios  c/D  =*  0.02 
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Table  la 

Commonly  Used  Thermocouples 


Types  of  Thermocouples 

Temperature 

Ranges 

°F 

°C 

Copper /Constantan 

-300/750 

-148/398 

Iron/ Cons tan tan 

-300/1600 

-148/871 

Chromel/ Alumel 

-300/2300 

-148/1260 

Chromel/Constantan 

32/1800 

0/982 

Platinum  10%  Rhodium/Platinum 

32/2800 

0/1537 

Platinum  13%  Rhodium/ Platinum  6%  Rh 

100/3270 

37/1798 

Platinel  1813  Platinel  1503 

32/2372 

0/1300 

Iridium/ Iridium  60%  Rhodium 

2552/3326 

1400/ 183(7 

Tungsten  3%  Rhenium/Tungsten  25%  Rhenium 

50/4000 

10/2204 

Tungsten/Tungsten  26%  Rhenium 

60/5072 

15/2800 

Tungsten  5%  Rhenium/ Tungsten  26%  Rhenium 

32/5000 

0/2760 

Table  lb 


Thermal  Properties  of  Thermalcouple  Materials 

cal  K3  *3  cal  P3C3  P3C3 

Metal,  Insulator  3 sec.  cm  < k steel  k aluminum  3 3 cmJ  pc  steel  pc  aluminum 
(subscript  3) 


Aluminum 

0.554 

4.66 

1.00 

0.65 

0.60 

1.0 

Copper 

0.935 

7.87 

1.69 

0.84 

0.78 

1.30 

Chromium 

0.201 

1.70 

0.36 

0.84 

0.78 

1.30 

Nickel 

0.160 

1.35 

0.29 

1.20 

1.12 

1.86 

Platinum 

0.174 

1.46 

0.31 

0.74 

0.69 

1.14 

Steel 

0.119 

1.0 

0.21 

1.08 

1.00 

1.66 

Tungsten 

0.373 

3.14 

0.67 

0.59 

0.55 

0.92 

Iridium 

0.345 

2.91 

0.62 

0.66 

0.61 

1.01 

Rhodium 

0.360 

3.03 

0.65 

0.72 

0.67 

1.12 

Rhenium 

0.111 

0.94 

0.20 

0.73 

0.68 

1.14 

Nickel-Chromium 

0.040 

0.34 

0.07 

1.01 

0.94 

1.56 

Copper-Nickel 

0.059 

0.50 

0.11 

0.99 

0.92 

1.52 

Teflon 

6 x lO-4 

0.005 

0.001 

0.54 

0.5 

0.83 

Air 

1.26  x 10~4 

0.001 

0.0002 

0.53 

0.49 

0.82 

1 

Note:  The  value  quoted  is  the  averaged 

value  over 

200  to  800 

°k  whenever 

data 

are  available. 
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PART  II  IMPROVED  ACCURACY  IN  THE  PREDICTION  OF  SURFACE  HEAT  FLUX 
AND  TEMPERATURE  BY  AN  INTRINSIC  THERMOCOUPLE 

1.  INTRODUCTION 

In  the  study  of  transient  heat  transfer,  many  experimental  diff icult  '.es 
may  arise  if  heat  flux  sensors  or  thermocouples  are  installed  direct  at  the 
surface  of  a body.  For  example,  a probe  may  be  damaged  by  a piston  or  a 
projectile  sliding  over  a cylinder  or  barrel.  A probe  on  a melting  and 
ablative  surface  of  heat  shield  can  be  easily  destroyed  because  of  high 
temperature.  Furthermore  a surface  probe  exposed  to  both  radiative  and 
convective  environment  may  measure  an  erroneous  surface  heat  flux  and 
temperature  if  the  probe  has  a different  radiative  property  from  that  of 
the  measured  surface.  In  these  circumstances,  calculation  of  the  transient 
surface  heat  flux  and  the  surface  temperature  can  be  achieved  by  inverting 
a temperature  history  measured  at  some  location  inside  the  body. 

In  general,  the  prediction  of  a surface  heat  flux  and  temperature  by 
the  measured  data  at  some  location  interior  to  a body  is  known  as  the 
"inverse  problem".  Many  configurations,  such  as  spheres,  cylinders,  and 
slabs,  had  been  studied  by  many  workers  and  many  methods  such  as  numerical, 
graphical,  series,  convolution  integral,  and  Laplace  transforms  were  used. 
Stolz  [1],  Beck  [2]  and  Williams  and  Curry  [3],  considered  the  numerical 
inversion  of  the  integral  solution  for  semi-infinite  and  other  bodies. 

In  this  method,  care  is  required  in  selecting  a time  interval  in  order 
to  achieve  a stable  solution.  Carslaw  and  Jaeger  [4],  Burggraf  [5], 
Koveryanov  [6],  and  Shumakov  [7],  respectively  considered  different  series 
approaches  in  which  generally  the  local  heat  flux  at  an  interior  location 
and  their  higher  derivaties  are  required.  However,  it  is  difficult  to 
measure  experimentally  or  to  process  the  measured  data  for  the  derivative 
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of  the  temperature.  Sparrow,  Haji-Sheikh,  and  Lundgren  [8],  Imber  and 
Kahn  [9],  Imber  [10],  Sabherwal  [11],  Masket  and  Vastano  [12],  Deverall  and 
Channapragada  [13]  and  Chen  and  Thomsen  [14]  applied  the  transform  method. 

In  these  works,  the  solution  is  represented  in  either  an  integral  form 
after  some  manipulation  of  the  contour  integral  from  the  inverse  trans- 
form, or  in  a series  form  after  an  expansion  of  the  solution  for  small  and 
large  times.  Using  Laplace  transformation  Chen  and  Thomsen  [14]  introduced 
a polynomial  in  terms  of  an  error  function  to  represent  the  response  of 
thermocouple  measurement  and  the  inversion  is  accomplished  for  any  transient 
surface  heat  flux  at  the  inner  surface  of  a cylindrical  tube.  In  thier 
study,  the  cylindrical  thickness  was  assumed  to  be  relatively  thick  such 
that  the  temperature  at  a large  distance  from  the  heating  surface  remains 
constant.  Therefore,  only  one  interior  temperature  response  near  the  sur- 
face was  needed  in  the  experimental  measurement.  Their  inversion  solution 
however  was  valid  only  for  a short  duration  due  to  the  asymptotic  expansion 
of  the  modified  Bessel  function  in  the  inverse  Laplace  transform.  Chen  and 
Chiou  [15]  studied  the  inversion  problem  for  the  case  of  a semi-infinite 
slab  or  a thick  slab  using  a Laplace  transformation.  The  exact  solution 
was  obtained  from  the  inverse  Laplace  transform  for  any  time  interval.  It 
was  then  shown  that  their  analysis  may  be  approximately  applied  to  the  case 
of  the  hollow  cylinder  if  the  interior  temperature  response  is  measured  at 
a location  close  to  the  inner  wall. 

This  report  presents  (a)  the  improved  numerical  solution  of  the  in- 
version solution  reported  by  Chen  and  Chiou  [15]  and  (b)  a further  demon- 
stration of  the  capability  of  the  solution. 
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The  theoretical  analysis  of  Chen  and  Chiou  [15]  is  recapitulated  in 
Appendix  I in  which  the  surface  heat  flux  and  temperature  is  predicted  by 
inverting  a temperature  history  measured  at  some  location  inside  the  solid 
body.  The  inversion  solution  is  obtained  by  invoking  Laplace  transformation. 
Both  the  surface  heat  flux  and  temperature  are  given  by  Eqs.  (19)  and  (20) 
in  Appendix  I. 

It  was  thought  that  the  accuracy  of  the  computer  program  generated 
for  the  solution  in  the  previous  report  by  Chen  and  Chiou  [15]  can  be  im- 
proved further  for  the  following  reasons.  First,  the  coefficients  b 

n 

(see  Appendix  Eq.  (11))  in  the  previous  formulation  has  a dimension  of 

temperature.  Therefore  the  determination  of  the  coefficients  depends  on 

the  temperature  range  of  each  particular  experiment.  It  was  found  that  the 

absolute  value  of  the  coefficients  b in  some  cases  can  become  as  large  as 

n 

44 

an  order  of  10  . Therefore  during  the  subsequently  numerical  manipulation 

in  the  computer  program  error  due  to  round  off  and  the  standard  fixed  up 
when  an  underflow  occurred  may  become  appreciable.  To  remedy  this  diffi- 
culty the  dimensionless  formulation  is  introduced  in  the  analysis  (Appendix 

I)  in  which  the  coefficient  b is  also  made  dimensionless.  As  a result 

n 

the  magnitude  of  the  coefficient  b can  be  greatly  reduced.  Secondly, 

n 

the  double  precision  format  was  not  used  throughout  the  previous  computer 
program.  It  is  felt  that  further  accurate  results  may  be  obtained  if  the 
double  precision  format  is  adopted  in  the  program. 

In  Section  II  the  new  computer  solution  is  shown  to  be  indeed  more 
accurate.  Also  in  Section  III  the  solution  is  shown  to  !be  capable  of  pre- 
dicting a case  involving  a periodic  surface  heat  flux  or  periodic  tempera- 


ture variation. 
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II.  RESULTS  OF  THE  IMPROVED  COMPUTER  PROGRAM 

The  previous  computer  program  of  Chen  and  Chiou  [15]  was  recasted  in 
dimensionless  form  and  written  in  the  double  precision  format.  The  new 
computer  program  is  listed  in  Appendix  II.  The  results  predicted  by  the 
new  and  previous  computer  program  are  given  in  Appendix  III  and  shown  in 
Figure  1 for  the  case  of  the  constant  surface  heat  flux.  This  is  the  case 
in  which  a steel  slab  initially  at  a uniform  temperature  is  suddenly  sub- 
jected to  a constant  heat  flux  Q at  one  of  the  surfaces  and  kept  at  the 
initial  temperature  on  the  other  surface.  Figure  1 shows  the  solution 
predicted  by  inverting  the  temperature  response  at  an  interior  of  the  slab 
from  the  new  and  previous  computer  program.  This  solution  predicted  by 
the  new  and  previous  programs  used  the  ten  term  representation  for  the 
thermocouple  response.  The  comparison  clearly  shows  the  improvement  of 
the  new  solution  over  the  previous  one.  Except  for  the  short  time  duration 
the  solution  with  the  new  program  reduces  the  error  to  only  one  half  of 
the  error  of  the  previous  program  i.e.,  an  error  of  less  than  one  percent. 

k 

In  the  short  time  period  the  solution  exhibits  a Gibbs  phenomenon  because 

I 

of  the  discontinuity  of  the  surface  temperature  gradient  occurred  at  initial 
condition.  The  solution  shows  a 17%  of  initial  overshot  of  heat  flux  and 
then  a 7.8%  of  undershot  before  the  solution  approaches  the  constant  heat 

flux.  It  should  be  remarked  that  Gibbs  phenomenan  is  artificially 

★ 

Gibbs  phenomena  [3]:  for  a sequence  of  transformation  Tn(t),  n = 1,  2,... 
of  a function  q(t)  (here  q(t)  = constant)  if  the  interval  * Tn(t), 

U“K)0 

lim  sup  Tn(t)]  contains  points  outside  the  interval  inf  q(t), 

t-*-t  n-»oo  t-^t 

lira0  ° 

t-*tQ  sup  q(t)]  then  the  sequence  is  said  to  exhibit  a Gibbs  phenomena. 
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introduced  due  to  the  idealization  of  the  initial  condition.  In  most 
practical  situations  the  surface  heat  flux  will  be  continuous.  Therefore, 

Gibbs  phenomenon  will  not  appear. 

Figure  2 (see  Table  1 also)  shows  the  comparison  between  the  solutions 
for  the  constant  heat  flux  case  with  10  and  20  term  representation  for  the 
thermocouple  response.  One  sees  that  the  solution  with  20  term  representa- 
tion after  the  initial  Gibbs  phenomenon  quickly  approaches  the  expected 
constant  heat  flux  solution  with  a negligible  error  of  less  than  0.14  per- 
cent. This  shows  the  accuracy  of  the  new  computer  program.  From  Figure  2 
one  also  observes  that  both  overshoot  and  undershoot  of  Gibbs  phenomenon 
are  smaller  for  the  20  term  representation.  Additionally  the  points  of  the 
overshoot  and  the  "ndershoot  has  moved  to  near  the  zero  time  which  agrees 
with  the  characteristic  of  Gibbs  phenomenon.  According  to  Gibbs  phenomenon 
the  point  of  overshoot  shoot  should  approach  the  initial  zero  if  the  number 
of  term  of  the  series  which  represent  the  thermocouple  response  is  increased 
to  infinite. 

VII.  VERIFICATION  OF  OSCILLATORY  SOLUTION 

As  a measure  of  applicability  of  the  present  inversion  solution,  a 
test  problem  was  solved  for  the  case  of  a slab  subjected  to  a preiodic 
surface  temperature  variation  on  one  surface  and  held  to  the  initial 
temperature  on  the  other.  The  analytic  solution  for  the  problem  is  given 
in  Appendix  IV  where  a more  suitable  form  of  the  solution  than  the  one 
given  by  Carslaw  and  Jaeger  [4]  is  derived  and  tabulated  for  the  thermo- 
couple response  at  one  tenth  of  the  slab  thickness  from  the  surface.  The 
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TABLE  1 

Comparision  of  Inversion  Prediction  and  Exact  Solution 


0.1  0.0624 

0.2  0.1628 
0.4  0.3395 

0.6  0.4882 

0.8  0.6183 

1.0  0.7353 


0(0, t) 

0(0,  t) 

ERROR'* 

80  (0,t) 

30  (0,t) 

ERROR 

EXACT 

PREDICTED 

-M. 

% 

3x 

EXACT 

8x 

PREDICTED 

% 

0.3758 

0.3773 

+0.399 

1.0000 

0.9873 

-1.27 

0.5315 

0.5314 

+0.019 

1.0000 

1.0014 

+0.14 

0.7516 

0.7516 

0.00 

1.0000 

1.0001 

+0.01 

0.9205 

0.9205 

0.00 

1.0000 

1.0000 

0.00 

1.0629 

1.0628 

-0.009 

1.0000 

0.9996 

-0.04 

1.1884 

1.1878 

-0.05 

1.0000 

1.9986 

-0.14 

0.3758 

0.4048 

+7.71 

1.0000 

1.0796 

+7.96 

0.5315 

0.5284 

-0.58 

1.0000 

0.9391 

-6.09 

0.7516 

0.7516 

0.00 

1.0000 

1.0179 

+1.79 

0.9205 

0.9205 

0.00 

1.0000 

0.9923 

-0.77 

1.0629 

1.0626 

-0.028 

1.0000 

1.0043 

+0.43 

1.1884 

1.1887 

+0.025 

1.0000 

0.9963 

-0.37 

TERROR  % = ((PREDICTED) -(EXACT)) /(EXACT) 
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surface  is  subjected  to  a periodic  temperature  variation  with  a period  of 
8 milliseceonds.  Fifteen  data  points  of  the  temperaturei  response  at  the 
thermocouple  location  are  then  input  to  the  inversion  program  for  prediction 
of  the  surface  temperature  and  heat  flux.  The  result  is  shown  in  Figure  3 
and  Appendix  III  where  the  data  symbol  "2"  denotes  the  thermocouple  re- 
sponse and  "1"  the  surface  temperature.  The  accuracy  of  the  inversion  pro- 
gram is  shown  in  Table  2.  Except  for  the  extremely  short  time  period  of 
0.4  milliseconds  the  prediction  by  the  inversion  program  with  15  term  repre- 
sentation is  within  2 percent  of  error. 

The  accuracy  can  be  improved  more  if  more  data  points  are  used. 

Figure  4 shows  the  predicted  surface  heat  flux  which  we  were  unable  to 
compute  from  the  series  solution  (Eq.  (5)  of  Appendix  IV).  This  demonstrates 
the  versatility  of  the  inversion  solution. 

IV.  APPLICATION  OF  THE  INVERSION  PROGRAM 

Three  sets  of  data  (see  Appendix  III)  provided  by  Rock  Island 
\rsenal  for  the  temperature  response  of  a thermocouple  embedded  in  a M60 
gun  barrel  were  utilized  to  evaluate  the  inversion  solution.  The  inversion 
prediction  for  the  surface  heat  flux  from  all  three  sets  of  data  were  ex- 
tremely high  when  compared  with  other  known  data  calculated  by  Chen  and 
Chiou  [15].  Since  the  program  correctly  predicted  the  surface  heat  flux 
for  other  sets  of  experimental  data  it  was  judged  that  the  three  sets  of 
data  may  contain  inaccurate  initial  time.  For  most  experimentations  the 
recording  instrument  is  likely  to  experience  some  delay  in  responding  to 
the  extremely  fast  transient  heat  flux  typical  in  gun  bores.  Therefore 
an  advanced  shift  of  time  of  2 millisecondc  in  the  data  was  tested.  The 
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Table  2 


Comparison 

of  Inversion  Prediction  and  Exact  Solution 

Time 

Surface 

Surface  Error A 

(sec. ) 

Temperature 

Temperature  % 

(theoretical) 

(predicted) 

0.0002 

189.5041 

193.4060 

+3.563 

0.0004 

296.3119 

301.4311 

+2.367 

0.0006 

397.7934 

403.8998 

+1.921 

0.0008 

491.4497 

498.3988 

+1.689 

0.0010 

574.9747 

582.5683 

+1.534 

0.0012 

646.3119 

654.3358 

+1.417 

0.0014 

703.7046 

711.9396 

+1.321 

0.0016 

745.7396 

753.9650 

+1.236 

0.0018 

771.3818 

779.3802 

+1.157 

0.0020 

780.0000 

787.5615 

+1.08 

0.0022 

771.3818 

778.3093 

+1.00 

0.0024 

745.7396 

751.8527 

+0.918 

0.0026 

703.7046 

408.8443 

+0.874 

0.0028 

646.3119 

650.3441 

+0.712 

0.0C30 

574.9747 

577.7933 

+0.569 

0.0032 

491.6697 

492.9790 

+0.372 

0.0034 

397.7934 

397.9893 

+0.0616 

0.0036 

296.3119 

295.2761 

-0.679 

0.0038 

189.5041 

188.8522 

-0.595 

*Error  = ((Predicted  - (Exact) )/ (Exact))  x 100 
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result  is  shown  in  Appendix  II1-3.  Figures  5,  6 and  7 respectively  show 

2 

the  predicted  maximum  heat  flux  of  1650,  1235  and  2695  Btu/ft  sec.  approx- 
imately at  2 milliseconds  after  the  firing,  while  Figures  8,  9 and  10  show 
that  the  surface  temperature  are  273,  234  and  396  degree  F.  approximately 
at  6 milliseconds  after  the  firing. 
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V.  CONCLUSION  AND  SUGGESTION 

The  new  Inversion  computer  program  was  thoroughly  tested  for  its 
applicability  and  accuracy.  The  program  can  invert  an  intrinsic  tem- 
perature response  to  predict  the  case  of  a constant  surface  heat  flux 
or  the  case  of  a periodic  surface  temperature  within  2%  of  deviation 
except  at  the  extremely  short  time  period.  This  was  achieved  with  the 
nm-ifiimtin  number  of  20  input  data  points. 

It  is  expected  that  the  accuracy  of  the  prediction  will  increase 
if  the  number  of  data  input  is  increased. 

Based  on  the  experience  gained  in  working  with  the  program  the 
following  suggestions  are  thought  relevant. 

(1)  In  conducting  an  experiment  it  is  vital  that  both  tempera- 
ture and  the  time  in  the  intrinsic  measurement  of  surface 
heat  flux  and  temperature  must  be  much  more  accurate  than 
the  direct  measurement.  This  is  because  the  inversion 
problem  always  involves  predicting  a large  heat  flux  or 
temperature  variation  from  the  data  with  small  variation. 
Therefore  a slight  error  in  the  time  or  temperature  measure- 
ment a large  error  will  result  in  the  prediction  of  the 
surface  heat  flux  or  temperature.  For  example,  in  the  pre- 
diction of  gun  barrel  heat  flux  considered  in  Section  III 
an  error  of  two  milliseconds  in  time  will  lead  to  100  per- 
cent error  in  prediction  of  the  surface  heat  flux. 
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(2)  In  selecting  the  data  points  for  input  to  the  computer  program 
care  must  be  excericised  not  to  create  a locally  abrupt  jump  in  the  data. 
An  abrupt  change  in  data  points  will  often  introduce  an  abnormal  fitting 
of  a curve  in  its  neighborhood  and  hence  resulting  in  an  incorrect  pre- 
diction of  the  surface  heat  flux  and  temperature.  If  indeed  the  abrupt 
jump  of  the  data  must  be  used,  then  more  data  points  in  its  neighborhood 
must  also  be  chosen. 
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APPENDIX  I.  ANALYSIS  OF  THE  INVERSION  PROBLEM 
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Consider  a slab,  having  a sufficient  wall  thickness,  L,  such  that  the 
outer  surface  temperature  has  a negligible  response  when  the  inner  surface 
is  exposed  to  a transient  heat  flux.  A probe,  for  example  a thermocouple, 
is  located  at  X =*  X^  and  it  is  normally  desirable,  as  reported  by  Chen  and 
Li  [16],  to  be  close  to  the  heating  surface  since  a better  transient 
response  and  more  accurate  experimental  measurements  can  be  obtained 
to  reduce  error  amplification  in  the  mathematical  inversion  program. 

Under  these  circumstances  we  thus  assume  in  the  analysis  L/X^  >>  1.  The 
governing  equation  for  the  transient  heat  conduction  may  be  written  in  a 
dimensionless  form  as 
2 

ae  a e n 

* 2 

3t  3x 

with  the  initial  and  boundary  conditions 


0(x,  0)  - 0 

(2) 

9(®,  t)  ■ 0 

(3) 

9(1,  t)  - f(t) 

(4) 

Where  x ® X/X^  is  a dimensionless  distance  from  the  heating  surface  and 
t = a t/X^  is  a dimensionless  time  or  Fourier  number  with  a being  the 
thermal  diffusivity  and  t the  real  time.  0 - (T  - Tq)/To  is  the  dimen- 
sionless temperature  above  the  initially  uniform  temperature  Tq  . 

f(t)  «*  (F(t)  - T )/T  is  the  dimensionless  measured  temperature  response 
o o 

at  x - 1 with  F(t)  being  the  measured  absolute  temperature.  The 


r 


I 


inversion  problem  is  then  given  the  interior  temperature  f(t)  to  predict 
the  surface  temperature  6(0,t)  and  the  surface  heat  flux  per  unit  area 
* -qW.OXj/T  it  . Here  k is  the  thermal  conductivity 

of  the  solid. 

The  above  problem  may  be  solved  by  Laplace  Transformation.  Let  the 
transformation  be: 


q(0,t)  or 


e(x,s) 


r 


9(x,t)e  ts  dt 


(5) 


where  0 is  continuous  otherwise  satisfied  the  Dirichlet's  condition.  The 
temperature  function  9 is  recovered  by  inversion  of  the  Laplace  Trans- 
formation as: 


9(x,t) 


•c+i° 


1 

2 rri 


9eSt  ds 


*c—  i® 


(6) 


where  c is  a suitable  positive  value.  Equation  (1)  and  the  initial 
condition  (2)  under  transformation  (5)  become: 


80 


(7) 


which  has  a solution,  with  integration  constants  A and  B, 


/s 


0(x,s)  • A e'SX  + Be 


-fs 


sx 


(8) 


The  transformation  of  boundary  conditions  (3)  and  (4)  into  the  Laplace 
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plane  give  0(®,s)  * 0,  0(1, a)  ■ f(s).  Substitution  of  the  boundary 
conditions  into  equation  (8),  we  get  0 and  its  derivative  as: 


0(x,s)  * f(s)  e 


i^(l  - x) 


- d “ ^ f(8)  e 


Ja(l  - x) 


According  to  Chen  and  Thomsen  [14],  we  choose  the  temperature  response  as 
measured  by  a probe  to  be  represented  by  a polynomial 


f(t)  - 51  b r(n+1)  i2"  erfc  (-M 

n=l  n 2A 


or  in  Laplace  plane 


f(s)  « e“^  £ r(n  + 1) 


The  b ' s are  coefficients  of  the  expansion  to  be  determined  such  that  the 
n 

N terra  polynomial  describes  the  temperature  response  f(t)  measured  at 
x • 1.  With  equation  (12) , equations  (9)  and  (10)  can  be  simplified  to 


0(x,s)  * e 


£r<"  + 1> 

n=l  s 


rx’/®  2 r(n  + 1) 


b 

n 

1/2+n 


The  inversion  of  equatior.s  (13)  and  (14)  at  x ■ 0 give: 


r 
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0(0, t) 


n 


n*l 


(15) 


N 


- ae(x,t) 


3x 


. V v n-1/2  r(tt4-l) 

x-0  SS  n r(n+l/2) 


(16) 


30 (0  t ) 

where  0(0, t)  gives  the  surface  temperature  and ' gives  the  surface 

qX  9x 

heat  flux  —r~  as  a function  of  time. 
icT 

o 


The  Integral  of  the  error  function  in  (11)  is  defined  as: 

I 

erfc  (y)  dy 


i^n  erfc(- 


(~)  - [ I2"'1 

2/7  -4-> 

2/t 


(17) 


f 

^ Jy 


-x 


with  n * 0,  erfc  (y)  ■ — 1 e dx  . r(n)  in  equations  (11)  and  (16) 

/v  J y 

is  the  gamma  function  or  Euler’s  integral  function  of  the  second  kind. 


T(n) 


r 

\ -ui  n-1 
1 e a) 

JO 


doj 


(18) 


It  should  be  remarked  that  the  choice  of  the  particular  form  (11)  is  to 

ensure  the  convergence  of  the  solution  on  the  Laplace  plane  and  an  analytic 

inversion  back  to  the  physical  plane.  With  b^  coefficients  determined 

from  equation  (11)  and  the  experimental  measurement  of  the  temperature 

response  f(t)  at  x * 1,  the  surface  temperature,  T (t),  is  obtained 

w 

from  equation  (15)  as 


N 

T (t)  - T (l  +22  b tnJ 
w ° \ n-1  n / 


(19) 


i 


3 HMD 
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-T  1C  30 

and  the  heat  flux,  q(t)  = ■ -°-  — 

dX 


x=0 


, from  equation  (16)  as: 


N 


, v -V  V . . n-l/2  r(n+l) 


(n+l/2) 


(20) 


The  above  solution  is  the  exact  solution  for  predicting  the  transient 
surface  heat  flux  and  temperature  valid  for  both  short  and  long  time  dura- 
tions as  long  as  the  slab  is  thick  enough  such  that  the  outer  surface 
maintains  its  initial  temperature.  The  feature  of  the  present  solution  is 
the  polynomial  (11)  which  on  the  Laplace  plane  gives  a term  in  equation 
(12)  exp  (-/£)  to  cancel  the  term  exp  (/s)  in  equations  (9)  and  (10). 
This  polynomial  (11)  as  suggested  by  Chen  and  Thomsen  [15]  makes  the  pres- 
ent solution  simpler  than  many  inversion  solutions  derived  in  the  past  and 
valid  in  any  time. 

If  the  fluid  temperature,  T (t),  away  from  the  surface  of  the  slab 

S 

is  known,  then  the  instantaneous  heat  transfer  coefficients,  h(t) , can 
be  determined  from  Newton's  cooling  law  as: 


h(t)  - 


3lL> 

T (t)  - T (t) 
g w 


(21) 


where  heat  flux,  q(t),  and  wall  temperature  T (t)  are  given  by 

w 

equations  (19)  and  (20),  and  the  average  heat  transfer  coefficient  up 

i 

to  time  t^  h(t)^  can  be  defined  as: 


h(t)  - 


C q(t')  dt’ 

Jo 


fC  [T  (t#)  - T (t ' ) ] dt’ 

Jn  8 w 


(22) 
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IMPROVED  COMPUTER  PROGRAM 


Cartesian  Inversion  Problem 


O O at  — « 
a«-Oiu c 
< « © & w • « 
W U O <1  O L 
O w X u/ 

x o ►-  o : 
w x • x z • v 
at  *-  © ►-  ui  o : 

Du.  x > bu  i 

U.  U.  z VO  • V 


lAOflf  CthKttlTWH 

cru«navi>oy* 


< UI 

ui  at  or 

UWD 
•20.1- 
• < X n 
*-  uj  a 
IV  K UI 

— — a 
o uj  r 

0 —J  uj 
>U1W  CL  H 

• at  -j  D 

- at  a o 'U 
. C “*occ 

/ (O  u c o 

1 b_J  X O 
■»  Sfi  o QL 

)3  t ivw 

• *-  a X < 

r © u ►-  o 
i < x 
r «iJ 


- »-  ~ z 

‘ r u<  z -• 

? O a — 

■ t/  N N 
I itfOO- 
r u»  m a — 
• *-v»  f »- 

) 1 - UIO 
3 V-  C H 


a*  o 

cr  «i 
a*  at  oc 
m at  b- 

iSb: 

O a/  ( 

• O'  30' 
O — LT  o O < 
NNHMN  ! 


«J  «l 

T V f3  f3  C I 
at  a ««  x «i  < 

CJ  C UJ  UJ  LU  I 

u.  u a at  at  i 


wOWOOOOOOOOOOOWWWOOOuOOUt 


O O O O O -u 


m « r^uiLP‘o-**'u*»0'*x>'*/r-«uvro-* 

rj  U 

DO  nO'J^O^OOO^DOOO') 

oo  000000000090000 


r 


0022  WR1TEI6.100J 

0021  100  FORMAT  I 20X,*  BORE  SURFACE  TEMPERATURE  A NO  FEAT  FlUX  PROGRAM* I 

00  2*  WRITE  It, 1091  MB 

002!  OC  10  1 = 1,  MP 

0024  REAOI 9 , 202  I TIME  1 1 1 .TEMP  1 1 ) 


NOW  e COEFFICIENT  FOR  THE  THE  TEMPERATURE  RESPONSE  AT 


“1  '•  O ♦ 

w . o o 

o ">  o o 

* • o 


a uIlU  J •<  w H 

Z O * ui  a r - 
* • ii  c j r o ui 

» o *n  ► ii  ii  w i r 

lOHJUi  — h N o 
'O-JUJ-Jf-O  — uj 
> • O 3 N ►. 

o o a — h 

h r*-  s:  M u>  *• 

i uj  uj  ■»  r* 

§C  M ar  o — 

a oiujk 
r> 


wia  o n <*»  / m 

•u  m ^ i«  i«  r» 

0*3000000 
OOWOOOOO 


■ONBV-O'FNW^Ifi 

oooooooooo 

OOOOOOOOOO 


« N HIT  O H N( 
lU  w W « CP  (M/> 

oocoooo 
oo  ouooo 


OHO  106  F OAMAT  I 1 NJKBER  OE  TIKE  TEMPERATURE  PAIRS  I SEC.  ,F.  1- • , I 21 

OKI  WRITE  C 6* 1 OS  I NR 

014]  WR1TEI6.107I 

014}  00  IS  1*1, NP 

0144  IS  WRITEI6.108I  TIKE  1 1 ) .TEMPI  1 1 


o 

o 


o 

o 


3 £ 


r 

Ul 

H- 

O 


II  *o 

o o 


cC  u>  ~D  O 

— ►- 

O - O < * 

►-  > -•  QC  • 

co  • mj  n 

• O a*  — 

a o n r «■* 

u»  u u.  uj  —• 

O h • H « 

o u.  k • 

— o * u. 

— I * x O 

►-  O m 

— mJ  m VI 

CD  -t  • - U* 

h-  O iu  ^ • 

ha  4 NI  UJ 

a O o a.  — •-  ■ 

X •►*•!»)• 

i/I  (X  U.  M • — — 


O K w 
Z <V 


K 

uj  s r> 


or  d o 
uj  a x 
a — — 

i5- 

X 


a — 


r-  » a.  es  O ►- 

* k u.  Z -J 

• P'S  ju  m 4 ^ 


; r -OOCOV 
/ u N O'  (M 

• Cw  — r*<  • 

H-  I -*  • * • • • 

• - « • n •»  *>  « 

— T — > 

• CL  ►-  ► * 

: i O < ••  < u uj  uj  i 

r u.  r «c  x h*  ►-  > 

a a tr  ~ ~ • 

> o o u a a it  a t 

u_  u.  On  I 3 2 ' 


• Z U 

X V3  QC 


• — O 
in  • uj 

• O Z 


S - 


■ O I — -»  *M  - «, 

1 r*>  a - • X 

* — z — — -o  - UJ 
i uj  r ►*  -- 


o — 
o • 

INI  'C 


► ■ u.  ►*  B I 


i < or  ou ; <i 

• r oub.i  o i 

• cf  - z "-cc  uj  r 
o -ictoauj' 

u • h J It  • H I 
•or -j  m 


-o  ► r — 

* ♦ ij  Ul  fl 

u.  UJ  k-  z 

W II  _ » 

M X H-  «0 

* uj  z w 

4 UJ  (l  U Ul 

r Hhoh 


4)  <J 

w w H 

UJ  UJ  o 


or  a O 
: j o 


t° 

it 


i- 


l x i 


o j sc 
O -j  X 
0^0 
-*  | 00 


J < <3 


a O uj 

X u.  O 


o 
o*  o 

O'  O 


O U o u)  o u> 


^I4)MU  JO**N4  JUI'O 

4 4 j I/I  «I  «I  •!  »i  »"T 


f»  9j  ir  o m r*  m -r  m 4MBVOHM 

m«n  IIHI  <J  41  4J  4J 


000000000000 


OO  OOOOOO  OOOOOOOO 


000 


0176  EN0 


CoS3  $UMTG=0.000 

OOSA  SUMO  =0.000 

00S5  SUMOT  = O.ODO 

CO't  TCI  11*1  TG(  1I-TEMPOI/(*60.000»TEHPO» 

C J * 7 00  80  1 = 1.200 


OOOOOOOOOOOOOOOOO 


WRITEI6.L0AI  ALP 

10A  FORMAT  <■  THERMAL  DIFFUSIVITT  IFTS0/SEC1  .'.FL5.81 
MAITEI6.105)  THCCN 

105  FORMAT!*  THERMAL  CCNC'JCTIVITTI  BTU/FT.SEC.F  . ) « ' ,F  15. 81 
WRITEI6.106)  NP 


f 


L- 

ft— 

< 

eg 

o 

• 

ec  'O 

u t u. 

i • 

r 

• x 

r* 

cc  ro 

• 

a • 

O 

—g 

uj  -o 

m4 

2? 

• • 

o 

ft 

X o 

ft 

Z 

a ^ 

— 

U.  UL 

X 

• 

ft 

• 

<* 

— 

II 

z 

UJ 

— X 

4 

o 

x trt 

a 

at 

VO 

w • 

o 

i- 

x in 

a 

» 

•J 

u.  • « 

CL 

*OUUJ 

X 

— *-<  a 

CL 

x u.  a z 

—ft 

w —» tu 

mi 

x m 

r 

VO 

■» 

** 

Z 

X 

X 

< 

ft— 

at 

— • o r x 

e- 

fft  i*  ** 

X 

»4*4  H 

y 

u 

«k  «■«  ft 

• CD 

z 

-ft  3 

ft* 

W —ft 

V0 

»- 

3 

UJ  4 UJ 

< cs  or 

O 

»-  X K 

t r?  i? 

a 

*4 

Mtf  M 

ec  co  »-  o 

® 

eg 

a o a 

O uj  uj  Z 

X U.  * 

U.OQCUJ 

VO 

-J 

UJ 

o 

> 

uh  r- 

ft* 

UJ 

«4  »4 

ft* 

> — o n 

l *>  • 

i eg  eo  t— 

I .o  - 

I -O  Ol  ft* 

I <M  w X 

I — k z x a 

i -4  ui  - • h 

II  ^ *4  H*S 

a a w ii  •«* 

I a ->  •-«  » 

I eo  uj  ~>  z 

I ♦ o O ct 
| D 

| < _J  < hO 

|uj  < O O a uj  Z 

| Ct  U O O ►*  C£  UJ 


I -J  I 'O 
r i eg 

1 u- 

v>  l eg 

UJ  | w ft 

2 |<I 

o ia^ 
1»~  1 


IC  I O fNI 

I a.  | eg  w 


N CO  W O <4 
m m «i  m m ^ f 

O O o “3  O O O 
o o o o o o o 


eg  eo  >r  tn  «o  r-  «J 
O «-»  o O O O W 

O O O o o o o 

O O O UJ  o o o 


pg  m ^ *r>  «o 

O u o g g 

o c*  o o o 

O O O <-J  <J 


00C7  00  l K= I* M 

C0C8  l OTQU  ,J)»QTQC  I,  JMTRANl!  ,K)*PHHK,  J » 

00CS  RETURN 

0010  END 


0001  FUNCTION  SPWUN,»,X,£P$,INOIC.NRC» 


62 


OO  V>w  w 
O O O O O 

uoooo 


o w*  o 

OUW- 

o o o o 
o vuo 


© o o 
o o o 


*-«-*-* 

© o © 
o o o 


**^>4NiyOlN««(VNN 

©OOOOOOOOOO 

ooooooooooo 


«D(A  O 
fM  oj  n« 
'JOT* 
O O O 


n r» 
O O 


O'* 

m m m m ^ # 

o o o o o o 
o o o o © o 


0JA2  17  IF(  J.NE.JC31RI  Alt,J»-A<I  , J )-A I JCK* A< I ROW*.  J1 

00*3  10  CCNTTNJE 

C ....ORDER  SClMTCN  VAUJf«  ANO  CREATE  J0R0  AARAT 


ECMAAM  IV  G LEVEL  21  SIMt'l  DATE  • 76157  15/31/00  PAGE  0002 


63 


APPENEIX  ITT  M6C  DATA 


M60  GUN  THERMOCOUPLE  7 15.0  INCHES  FROM  BREECH 
BURE  SURFACE  TEMPERATURE  AND  HEAT  FLUX  PROGRAM 
NUMBER  UF  B(IJ  CCFF.  TO  BE  FITTED  * 10 
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APPENDIX  IV.  THE  CASE  OF  OSCILATORY  SURFACE  TEMPERATURE 


Consider  a slab  with  a sufficient  thickness,  i,  such  that  when  a surface 
is  subjected  to  a periodic  surface  temperature  variation  with  a frequency 
w the  other  surface  is  held  at  the  initial  temperature  T . If  properties 
are  assumed  constant  the  governing  equation  for  the  problem  can  be  written 


with  initial  and  boundary  conditions  as 


v(x,  o)  « 0 


v(£,  t)  ■ sin  wt 
v (o , t)  * 0 


where  v ■ (T  - T )/(T  - T ) and  a is  the  thermal  diffusivity. 

o max  o 

The  solution  of  the  problem  according  to  Carslaw  and  Jaeger  [A]  can 
be  written  as 


oo  r 

■ 2afT  ^ (-1)  (-l)nn  L(an  tt2  sin 


2 2 ,„2. 


wt  - w£  cos  wt) 


. /niTx\  ,.244,  2 2, 

• sin  (“yj  /[a  n ir  + w f,  ] 


Part  ITT  PREDICTION  OF  TRANSIENT  SURFACE  HEAT  FLUX  AND 
TEMPERATURE  ON  A HOLLOW  CYLINDER 

I.  INTRODUCTION 

In  the  study  of  transient  heat  transfer  many  efforts  have  been  made 
on  the  so-called  "inverse  problem"  [1,2]  where  a surface  heat  flux  and 
temperature  is  to  be  predicted  by  the  measured  data  at  some  location  in- 
terior to  a body. 

In  the  previous  works  [1-6]  the  solution  is  represented  in  either  an 
integral  form  after  some  manipulation  of  the  contour  integral  from  the 
inverse  transform,  or  in  a series  form  after  the  expansion  of  the  solution 
for  small  and  large  times.  Using  Laplace  transformation  Chen  and  Thomsen 
[6]  introduced  a polynomial  in  terms  of  the  error  function  to  represent  the 
response  of  thermocouple  measurement  and  the  inversion  is  accomplished 
for  any  transient  surface  heat  flux  at  the  inner  surface  of  a cylindrical 
tube.  However,  their  inversion  solution  is  valid  only  for  a short  duration 
due  to  the  asymptotic  expansion  of  the  modified  Bessel  function  in  the  in- 
verse Laplace  transform.  In  this  study  an  exact  solution  obtained  from  the 
inverse  Laplace  transform  by  the  convolution  method  is  given  for  the  case 
of  hollow  cylinder.  The  solution  is  valid  for  both  constant  and  variable 
heat  flux  and  for  both  short  and  long  time  duration. 

II.  ANALYSIS 

Consider  a long  hollow  cylinder  with  sufficient  wall  thickness  such 
that  the  outer  surface  temperature  has  a negligible  response  when  the  inner 
surface  is  exposed  to  a thermal  pulse  of  a transient  process.  This  condi- 
tion considerably  simplifies  the  theoretical  analysis  as  the  outer  boundary 
may  be  assumed  to  be  infinite,  and  only  one  interior  probe  of  the  cylinder 
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is  required  in  the  experimental  measurement.  The  material  of  the  cylinder 

» 

is  considered  to  be  homogeneous  and  isotropic  with  constant  thermal  dif- 
fusivity,  o.  Let  and  Rq  be,  respectively,  the  inner  and  outer  surface 
radii.  the  radius  of  the  probe  location  and  t the  dimensionless  time. 
If  the  temperature  of  the  cylinder  is  initially  uniform  at  Tq,  the  mathe- 
matical problem  governing  the  temperature  T,  may  be  written  as 


30  _ 3^0  1 38 

3t  " 9r2  r 9r 


1 < r < r * 
o 


(1) 


9 (r,  o)  - 0 

9 (o»,  t)  * 0 


(2) 

(3) 


0 (r^,  t)  ■ f(t)  1 < r^  < ® 


(4) 


7 

where  0 ■ T - Tq,  r = R/R^,  t * ccr/R^  , and  f(t)  is  the  interior  tempera- 
ture response  of  the  thermocouple  measured  at  r * r^  at  the  dimensionless 
time  t.  The  problem  is  to  predict  the  surface  temperature  0 (1,  t)  and 
heat  flux  per  unit  area 


q - - (K/R  )(39/3r) I 

'r  = 1 


(5) 


where  K is  the  thermal  conductivity. 

The  problem  can  be  solved  by  Laplace  transformation.  Let  the  trans- 
formation be 

0e“tSdt  (6) 


9(r,  s)  - f 


I 


when  6 satisfies  the  Dirichlet's  condition  the  temperature  function  9 is 
recovered  by  inversion  of  the  Laplace  transformation  as 
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8<r’  " 2^1 


-c+i«° 

•c-i* 


9eStds 


(7) 


where  c is  a suitable  positive  value.  Equation  (1)  and  (2)  under  trans- 
formation (6)  becomes 


d2e  . l de  - 

— - + -—-=*  s9 
.2  r dr 
dr 


which  has  a solution  of  the  form 


0 « AI  (pr)  + BK  (pr)  (8) 

o o 

where  I and  K are  modified  Bessel  functions  of  the  first  and  second  kind 
o o 

1/2 

with  p = (s)  . With  the  boundary  conditions  (3)  and  (4),  (8)  becomes 

6 - f(s)  [Ko(pr)/Ko(pri)]  (9) 

where  f(s)  is  the  Laplace  transform  of  the  boundary  condition  (4). 

The  temperature  response  measured  at  r * r^  can  be  expresses  by  a 
polynomial  or  numerous  other  suitable  functions.  In  the  present  analysis, 
for  reasons  to  be  explained  later,  f(t)  will  be  represented  as 


f(t) 


N -t 

£ ^ / ’I 

n-1  *o  „ 


(t)  Fn(t  - x)dr 


If  we  choose  F^(t)  * — e 


11 

4t 


(10) 


and  Fn(t  - t)  being  any  arbitrary  function 


I 


depending  on  n,  for  example  (t  - r)n  etc.  then  the  Laplace  transform  of 
Eq.  (10)  gives 


f(s)  - £ bn  Ko(pri)Pn(s) 


where  Fn(s)  is  the  Laplace  transform  of  Fn(t).  Substituting  Eq.  (11)  into 
Eq.  (19)  we  have 


5(s,  r)  - ]T  bn  KQ(pr)Fn(s) 


It  is  noted  that  the  K^(pr^)  in  Eq.  (9)  has  been  cancelled  by  this  sub- 
stitution which  explains  the  choice  of  F (c)  = ^ e__^  in  Eq.  (10). 
The  inversion  of  Eq.  (12)  gives 


N „t 


0(t,  r)  - £ bn  f j-e  4T  Fn(t  - T)dT 
n=*l 


At  surface  r ■*  1 


9 (t , 1) 


N rt 

R / £ 


e Fn(t  - x)dt 


The  temperature  gradient  and  hence  heat  flux  at  surface  is 


“ b f e 4t 

r-1  n *4  4t2  Fn(t  • T) 


some  examples  of  Fn(t  - t)  frunction  and  the  representation  of  the  thermo- 
couple response  are: 


I 
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Case  1 


Case  2 


Case  3 


Case  4 


Case  5 


If  Fn(s) 


2 2 
s +n 


N -t 

then  f(t)  - L bn  J Y7 
n»l  •'o 


ri 

4t  1 


— sin  n(t 
n 


- t ) dx  (16) 


If  Fn(s) 


2 2 
s +n 


then  f(t) 


N ft 

£ bn  l TT 


1 

4t 


cos  n(t  - x)dx  (17) 


If  Fn(s) 


then  f(t) 


(s  + a)1 


N rt 

n=l  */o 


~ 4x  >n-l  -a(t-x) 

e ct2;(7-  if — (18) 


If  Fn(s)  * s 


-(n+1/2) 


N -t  - 


t(t) . f;  f .-*■ 

n=l  •'o  2x  1.3.5. . . (2n  - 1)/tt 


then 


dx  (19) 


If  Fn(s) 


then  f(t) 


N rt 

/ 

n*l  •'o 


_ll_ 

4x  . .n-1 

e (t  ~ x) 


2T(n  - 1)  l 


r dx 


(20) 
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